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The continuum Yang-Mills functional integral within the first order formalism and in Coulomb 
gauge is studied. In particular, the temporal zero-modes of the Faddeev-Popov operator are explic- 
itly accounted for. It is shown that the treatment of these zero-modes results in the constraint that 
the total color charge of the system vanishes at all times. Further, it is argued that the functional 
integral is effectively fully gauge-fixed once Gauss' law has been resolved in Coulomb gauge. 
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O . 1. INTRODUCTION 

<n ; 

Coulomb gauge quantum chromodynamics [QCD] and Yang-Mills theory have recently begun to attract consid- 
erable attention, primarily because of their great potential in studying confinement. This potential has long been 
recognized and led to the Gribov-Zwanziger scenario for confinement [1-3] whereby the temporal component of the 
gluon propagator provides for a long-range confining force whilst the transverse spatial components are suppressed 
in the infrared. However, progress in Coulomb gauge has been hindered by the inherent noncovariance of the gauge 
condition. Gratifyingly, the technical obstacles are being steadily overcome. 

Various approaches to Coulomb gauge are currently being considered. There exist lattice studies [4-9] that are 
beginning to shed light on the nonperturbative behavior of the propagators. Also, a Hamiltonian-based approach to 
the problem [10-13] based on the original work of Ref. [14] describes various features of the system in a coherent 
J^" 1 , way [15-17]. The Lagrange-based (Dyson-Schwinger) approach to the problem [3, 18] is also making progress. One 
common theme to all these studies (and central to the Gribov-Zwanziger confinement scenario) is the importance of 
i—i, the Faddeev-Popov operator that arises after one fixes the gauge. 

In addition to its noncovariant nature, Coulomb gauge is also incomplete — the gauge is only partially fixed. After 
applying Coulomb gauge, one can still perform time-dependent (spatially independent) gauge transformations and this 
raises potential questions about the strict validity of Coulomb gauge [19]. On the other hand, trying to completely fix 
£f*l I the gauge in the continuum seems to lead to a contradiction [20] . Related to the gauge-fixing and noncovariance issues 
. in Coulomb gauge is the problem of energy divergences. Because the gauge-fixing only involves spatial derivatives, 
the ghost propagator is independent of the energy and closed loops involving these ghosts appear at first sight to 
exhibit a pure (and unregularizable) energy divergence. It has been argued that these divergences cancel to all orders 
perturbatively [21]. This cancellation, which has been explicitly verified at one-loop order [22, 23], is not arbitrary: 
within the first order formalism, the Faddeev-Popov determinant cancels when the temporal gauge field is integrated 
O ■ out [3, 18]. 

In this paper, the Faddeev-Popov operator is considered with the focus on the temporal zero-modes. These modes 
embody both the noncovariant and incomplete aspects of the Coulomb gauge- fixing. The first order, (Lagrange-based) 
functional formalism in the continuum is employed. It is shown that after Gauss' law is implemented (which in the 
functional approach arises after integrating out the temporal gauge field), the temporal zero-modes can be accounted 
for and their resolution leads to a completely gauge-fixed action that is free of energy divergences (in the sense that 
the Faddeev-Popov determinant cancels) with a well-defined, conserved and vanishing total color charge. 

2. FUNCTIONAL FORMALISM 

The starting point for this study is the functional integral associated with continuum Yang-Mills theory (see Ref. [18] 
for a complete description), 
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Z = Jv® cxp {iS YM }, Sym = ~ J d 4 x [E 2 - B 2 ] 



(2.1) 



with the chromoelectric and chromomagnetic fields given in terms of the spatial (A) and temporal (a) gauge fields: 

E a = -d°l a -Va a + gf abc A % a c = -d°A a -D ac [A}a c , 

B? = e ijk ~<?r bc 4A|] (2.2) 
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and where 

Df[A] = 6 ab V x - gr b A% (2.3) 

is the spatial covariant derivative. 2?<I> denotes the functional integral measure for all field types that may be present. 
The action is invariant under local (and global), infinitesimal gauge transformations of the form: 

Af = A% + - g D?[A]e b x , 

of = <-- g dX-f acb * L X- (2-4) 

Notice that the functional integral measure, 2?$, includes all those configurations related by such gauge transforma- 
tions and thus contains the integration over the gauge group. Since the action is gauge invariant, integrating over 
the gauge group produces a (divergent) global factor which, in principle, can be absorbed into the normalization of 
the functional integral and in itself is harmless. However, inconsistent integration over the gauge group does become 
problematic when calculating gauge-variant Green's functions from the functional integral: the integration over the 
gauge group averages such quantities to zero. For example, in the free (non-interacting) theory it is necessary to fix 
the gauge in order to be able to invert the quadratic part of the action such that particle propagation may be properly 
defined. We fix the gauge by applying the Faddeev-Popov technique. This involves inserting the following identity 
(9 X is the parameter associated with the gauge transformation, see Eq. (2.4)) into the above functional integral: 



1 = J m x S(F[A])det [M ba (y, x)] , M ba (y, x) = 



(2.5) 

F=0 



We are concerned here with Coulomb gauge: 

F a [A] = V-A a = (2.6) 

and for which the Faddeev-Popov kernel reads 

M ba (y,x) ~ -V x -Df[A]5(y-x). (2.7) 

There are caveats to the identity, Eq. (2.5), above. When the gauge fixing is incomplete, zero-modes of the Faddeev- 
Popov operator will arise, i.e., 9 or A are such that F[A e ] = F[A] and for Coulomb gauge, Eq. (2.6), there are two 
cases: 

1. — Vx-D^IA^ — 0, Gribov copies for A ^ (generated by spatially and temporally dependent gauge transfor- 
mations) 

2. — V x ■ D x b [A]9 b {t) = 0, temporal zero- modes (generated by temporally dependent but spatially independent 
gauge transformations). 

[We implicitly include global (6* = constant) transformations into the latter category.] Clearly, in both these cases, 
M has zero eigenvalues and consequently the Faddeev-Popov determinant also vanishes, violating Eq. (2.5) since the 
left-hand side cannot vanish. We must therefore modify the formalism to account for this and we write 

1 = J V9 x S(F[A])drt [M ba (y,x)] , (2.8) 

where T>9 X explicitly excludes any spatially independent 9(t) and 

del [M ba (y, x)] = dct [M ba {y, x)] (2.9) 

is the determinant with the zero-modes of the operator (temporal or Gribov copy) removed. We can now write our 
functional integral as 

Z = J V^S(V-A a )det -V-D cxp{iS YM } (2.10) 

where it is understood that the functional integration measure, £>$, still contains the integration over the full gauge 
group. The direction(s) in group space corresponding to the zero-modes of the Faddeev-Popov kernel are still explicitly 
present within the functional integration and may still cause problems. 
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To proceed, we convert to the first order formalism [3, 18]. We introduce an auxiliary vector field (7?) via the 
following identity: 



cxpji J d 4 x^E a -E a ^ = J 2?7rexp|z J d 4 x 



7f a -7f a -TT a -E a 

2 



The 7T-field is further split into components using the identity 



const = J V(j)VT expj-7 J d 4 xr a (V-7r a + V 2 a ) J, 
changing variables 7? = 7? — V<fr and integrating out the Lagrange multiplier field, t, to arrive at the form 



/ 



2?$<5(V-A a )(5(V-7T a )dct 



-V-D 



exp {iS} 



where 



J d 4 x -±B 2 -^n 2 -^(V<t>f+TT a -d°A a + a a {v-D^cf + gp a ) 



and with the color charge of the gauge field 



pa = jade A d_ r 



(2-11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



The advantage of the first order formalism is, of course, that the action is linear in a and this field can be integrated 
out to give 



Z = J 2?$(5(V-i a )(5(V-7f a )det S (v-D ab <t> b + .gp a ) exp {iS'} 



where 



^ r>2 1 _2 1 en i\2 



IB* - ^n z - ^{V(f>Y + n a -d°A a 



(2.16) 



(2.17) 



The ^-functional constraint on the scalar field (j) is the functional expression of Gauss' law. To resolve the constraint 
and eliminate the functional integration over the scalar field (f>, we first must investigate the Faddeev-Popov operator 
a little more closely. The eigenvalue equation for the operator is 



-V-Df[l]^(f;t) = A n <(f;i) 
with a complete orthonormal basis <fin(x; t) satisfying 

d 4 xip™{x;t)ip a n {x,t) = S mn . 



(2.18) 



(2.19) 



However, since the time dependence of the operator is only implicit (within the spacetime dependent field A x ), the 
time argument is only a label (in the sense that there is a different spatial operator at each time). We could state that 
the eigenvalue equation is evaluated at a specific time and for which the eigenfunctions are spatially orthonormalized 
at this time. In this case, the eigenvalues become implicitly time dependent in that they refer to the eigenvalue 
equation at a particular time, t: 



-V-Df[A} v b tn (x) = \ tnV a t 



a tn{x), jdf 



XifZixMnix) = 6 V 



(2.20) 



We use the convention that \tn=o = denotes the collection of zero modes of the Faddeev-Popov operator. The 
temporal zero-modes are spatially constant and for SU(N C ), there are N 2 — 1 such linearly independent eigenvectors 
in the color space which when necessary, we label by an additional index fi: (fito(x) = ip b ^ (the adjoint color index, b 
here, labels the component of the /i-th vector). Using the complete orthonormal basis, Eq. (2.20), we can expand the 
field, 4>, and color charge, p, as series (index /1 is implicit within the zero modes): 



(2.21) 



71 = 



n=0 
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where the coefficients are time dependent and, in particular, 

ao„(t) ~ J d 3 Xip* t ^{x)p a x . 



(2.22) 



Having identified the zero-modes, albeit formally, we may define the inverse Faddeev-Popov operator in their absence, 
M: 



-V-D 



M* T = M 



-V-D 



such that 



detM = Y[ Kn = det 



-V-D 



(2.23) 



(2.24) 



Returning to the functional integral over (f> within the expression, Eq. (2.16), we may now write 

1 



J V<f>5 (y-D ah 4> b + gp a ^j cxp 1 1 J d A x 



Y[ db n (t)S (-b n (t)\ tn + ga n {t)) 



-My 



.71 = 



m,n— 



Y[S(ga ^t)) [ dbo^t) [ 

., J — OO J 



jj db n (t)S {-b n {t)X tn + ga n (t)) 



exp 



|6 0M (t)| 2 ^(f)V^S, M (£) + 2 6S M (t)6„(t)^(f)V^? n (£) 



n=l 



£ C(t)6„(t)^(f)V^?„(f) 



[j5(5ao M (i))det 



-V-D 



d6 0At (i) x 



exp • 



d 4 x 



|6 ,(t)| 2 ^(f)V^S, M (f) + 2 ff f; &5 M (*)¥^(*)V^(f) 



+5 x * x ^tmWVr^tnW 



m,n=l Lnt 



(2.25) 



If we neglect the (spatially dependent) Gribov copies, the temporal zero modes can be completely eliminated by noting 
that ^xftOfi = 0, such that the integrals over the bon(t) are overall (divergent) constants which can be absorbed into 
the normalization [26]. Thus, using Eq. (2.22) 



J V4>5 (V-D<V + 9P a ) exp |t J 



d 4 x 



~<y*r 



-V-D 



exp 



d 4 x 



9 2 £ a h { P a : {t K zm 2 M^ 



— j A t * m A t „ 



. (2.26) 



Now, because the temporal zero- mode is also a zero-mode of the Laplacian, we can infer that the ratio V^to^/Ato 
is finite and that if we multiply by a 0fl (t) — then the product vanishes. Thus, the sum in the exponential can be 
extended to include the zero-modes without reintroducing any ambiguities (we are, after all, only adding zero) and 
we have the form 



J v<t>5 (i 



V-D a V + 5P a exp^i / d*x 



d 3 xp° c I det 



-V-D 



exp { - I d A x 



-V X -D?[A] 



-V X -D*[A] 



(2.27) 
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where the quantities in the exponent formally include the zero-modes. As demonstrated, these zero-modes play no 
role other than to allow us to write the action conventionally in terms of the full fields; however, nontrivially, these 
full fields have been shown to be well-defined insofar as the temporal zero-modes are concerned. Since the ^ 0fl are 
(Af? — 1) linearly independent vectors in the adjoint color space, we also have that 



(2.28) 



Returning to the original functional integral, Eq. (2.16), and substituting in the above result for the <j> integration, 
Eq. (2.27), we see immediately that the modified Faddeev-Popov determinants cancel, leaving 



W-A a )S(V-Tf a )d d 3 xp a x )exp{iS"} 



with 



S" = 



d A x 



--B 2 - -7r 2 + 7r a -d°A a + -g 2 p\ [-V X -D?[A] 



Px 



(2.29) 



(2.30) 



Given the definition of the color charge p a , Eq. (2.15), it is clear that the argument of its 6- functional constraint 
(J d 3 'x f abc A 6 ■ tt c ) must vanish at each time, t. Further, the vanishing of the charge is invariant under temporal 
gauge transformations since (using the fundamental representation for the colored fields and denoting the spatially 
independent gauge transformation in this representation U t ) 







J d 3 xp = J d 3 x A,n -> J d 3 x U t AU},U t nU} =U t J 



A,TT 



U 



0. 



(2.31) 



The action, S" given by Eq. (2.30), is however no longer invariant under such temporal transformations (courtesy of 
the TT a -d°A a term which remains after we have integrated out the tr-field) and the temporal zero- modes of our new 
functional integral, Eq. (2.29), do not exist. This functional integral is thus fully gauge-fixed (it contains no more 
zero-modes, except those Gribov copies that we ignore here) and constrains the conserved, total color charge to be 
vanishing: 



J d 3 xp a x = J d 3 xf abc A b -7r c 



0. 



(2.32) 



The appearance of the ^-functional constraint in Eq. (2.29), which ensures the vanishing of the total color charge is 
the main result of this work. Given that the action, Eq. (2.30), is unchanged from the form originally derived in [3], 
the consideration of the zero-modes leads to the further observation that the dynamics of the theory are unaltered. 
If quarks were to be included in the original action, then the color charge acquires an extra component: 



p a q =qT a !°q 



(2.33) 



where T a is the (Hermitian) generator of the gauge group, q is the quark field and q its conjugate. The total charge 
is constrained as before, i.e., 



J d 3 x [ P a + p a q ] =0. 

The inclusion of quarks into the first order formalism is currently being studied [24] . 



(2.34) 



3. SUMMARY AND CONCLUSIONS 



The temporal (and global) zero-modes of the Faddeev-Popov operator for Coulomb gauge, continuum Yang-Mills 
theory within the first order formalism have been studied. The explicit separation of such modes leads to a well-defined 
total color charge that is conserved and vanishing at all times. 

The appearance of the total color charge is interesting for several reasons. Firstly, that the total charge of the 
system is conserved and vanishing is a necessary condition for confinement — one could not expect only color singlet 
hadrons to emerge if the system were not so. As such, the properties of this charge may be regarded as a proof that 
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at the least, the total system (i.e., the universe) is colorless. This does not however say much about the confinement 
of quarks and gluons, or the observed spectrum of hadrons. That being said though, whilst the temporal zero-modes 
lead only to a total conserved color charge, one may speculate about the role of the Gribov copies (i.e., the spatially 
dependent zero-modes) in confinement. This is the premise of the Gribov-Zwanziger confinement scenario [1-3] and 
the above analysis intuitively supports this. 

Second, given that the temporal zero-modes drop out of the functional integral without modification of the effective, 
gauge-fixed action despite the fact that the (Coulomb) gauge-fixing is incomplete leads to the conclusion that no 
further gauge-fixing considerations are necessary and that Coulomb gauge is a well-defined choice of gauge insofar as 
the dynamics are concerned (leaving aside the definition of the physical state space). Further, since such temporal 
zero-modes are associated with the ghost energy divergence problem it is clear that these divergences must cancel in 
final (physically meaningful) expressions, at least in principle (although they may occur in individual components). 

Thirdly, Gauss' law (which here appears as the ^-functional constraint on the auxiliary </>-field in Eq. (2.16) and 
arises when one integrates out the temporal, cr-field) plays a pivotal role in the extraction of the zero-modes and 
therefore in both the definition (and conservation) of the total color charge and in resolving the incompleteness of the 
gauge (see the previous paragraph). This is not surprising, since Gauss' law defines the charge via the generator of 
gauge transformations [16]. In the Hamiltonian approach to Yang-Mills theory [14], the imposition of Gauss' law as 
a constraint on the physical state space ensures gauge invariance in exactly this way. The resolution of Gauss' law in 
Coulomb gauge results in the explicit appearance of the nonabelian Coulomb potential in the Hamiltonian (and also 
here as the last term in the effective action, Eq. (2.30)). 

Fourthly, whilst the first order formalism has been employed here, the manipulations required to construct this 
formalism from the more usual second order formalism comprise identities that merely serve to linearize the action 
with respect to the temporal gauge field, a. Thus, what is true in the first order formalism will also be true in 
the second order formalism, although it will be manifested differently. The eventual form of the functional integral, 
Eq. (2.29), with the effective action, Eq. (2.30), here are non-local expressions and as such are not amenable to 
renormalization (that the first order formalism is not multiplicatively rcnormalizable is known [3]). In the standard, 
local version of the second order formalism (i.e., where the Faddeev-Popov determinant is expressed in terms of 
Grassmannian ghost fields etc.) which appears to be locally renormalizable (although there is as yet no conclusive 
proof) and therefore more suited to calculation, it may indeed not be possible to write the constraint that the total 
color charge is conserved and vanishing in convenient form but this does not mean that this is not true — quite the 
opposite. In any event, the dynamics of the theory are unaltered and in particular from a practical standpoint, the 
incomplete gauge-fixing poses no fundamental problem and the ghost energy divergences will cancel in the second 
order formalism just as in the first. 

It is worth pointing out that in principle, the 7?-field in the final expressions, Eq. (2.29) and Eq. (2.30), can be 
integrated out (the exponent is at most quadratic in 7?) and the ^-functional constraint on the total charge explicitly 
resolved. Whether or not the resulting expressions have any practical use however, remains to be seen, because of 
their highly non-local nature. 

Finally, the gauge invariance of the theory as applied to the Green's functions leads to Ward-Takahashi and Slavnov- 
Taylor identities. If one naively applies a temporal gauge transformation and tries to construct these identities, the 
resulting expressions are meaningless since the functional integration measure is not properly defined. As such, this 
is the only visible effect of leaving the temporal zero-modes of the Faddeev-Popov operator in the functional integral. 
However, once recognized this poses no problem in practice since to derive such identities one generally considers a 
fully spacetime dependent gauge transformation. The derivation of such identities is the focus of present work [25] . 
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